The Abelian Manna model on two fractal lattices 
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We analyze the avalanche size distribution of the Abelian Manna model on two different fractal 
lattices with the same dimension d g = In 3/ In 2, with the aim to probe for scaling behavior and 
to study the systematic dependence of the critical exponents on the dimension and structure of 
the lattices. We show that the scaling law D(2 — r) = d w generalizes the corresponding scaling 
law on regular lattices, in particular hypercubes, where d w = 2. Furthermore, we observe that 
the lattice dimension d g , the fractal dimension of the random walk on the lattice d w , and the 
critical exponent D, form a plane in 3D parameter space, i.e. they obey the linear relationship 
D = 0.632(3)d 9 + 0.98(l)d™ - 0.49(3). 



PACS numbers: 05.65.+b, 05.70. Jk, 64.60.Ak 

Although extensive research has been performed on 
self-organized criticality [1] for models on hypercubic lat- 
tices, far less work has been done on fractal lattices [3] . 
It remains somewhat unclear what to conclude from the 
latter studies. Fractal lattices are important for the un- 
derstanding of critical phenomena for a number of rea- 
sons. Firstly, results for critical exponents in lattices with 
non-integer dimensions might provide a means to deter- 
mine the terms of their e = 4 — d expansion. Secondly, 
fractal lattices are particularly suitable for a real space 
renormalization group procedures, in particular that by 
Migdal and Kadanoff [IHSj- Thirdly, scaling relations 
that are derived in a straightforward fashion on hypercu- 
bic lattices can be put to test in a more general setting. 
In this Brief Report, we address the first and the third 
aspect, by examining both numerically and analytically 
the scaling behaviour of the Abelian version of the Manna 
model [THHj on two different fractal lattices. 

The fractal lattices used in this study are generated 
from the arc-fractal system [10 . The lattice sites are 
the invariant set of points of the arc-fractal. We con- 
sider nearest neighbor interactions among sites. Here, 
the nearest neighbors of a given site are all sites which 
have the (same) shortest Euclidean distance to it. Our 
fractal lattices have no natural boundary; instead, they 
have only two end points at which two copies can join to 
form a bigger lattice. The dimension of the lattices is the 
same as the arc-fractal that generates them. 

In this study, we shall consider two fractal lattices: 
the Sierpinski arrowhead and the crab (see Fig. [I]). The 
former is named "Sierpinski arrowhead" because it is 
the same as the well-known Sierpinski arrowhead |llj . 
whereas the latter is termed "crab" because the overall 
shape of the generated lattice looks like a crab. These 
fractal lattices are generated through the arc-fractal sys- 
tem with number of segments n = 3 and opening angle 
of the arc a = n. For the Sierpinski arrowhead, the rule 



for orientating the arc at each iteration is "in-out-in" , 
while the rule is "out-in-out" for the crab. Both lattices 
have the same dimension d g = In 3/ In 2 ss 1.58. The to- 
tal number of sites on the lattice at the i-th iteration is 
Ni = 3 l + 1. The coordination number of sites on these 
lattices varies between two and three. Asymptotically, 
one-third of the sites have three nearest neighbors (called 
extended sites) , while the remaining two thirds have two 
nearest neighbors (called normal sites). Since the lattice 
sites are being stringed up by arcs (see Fig. [I]), they can 
be labeled as sites on a ID linear chain. In order to de- 
termine the linear size (to be used in finite size scaling) 
of the lattice, reference sites (hollow circles, which are 
not part of the lattice) have been added between the real 
sites (solid circles). This is possible because of the uni- 
form spacing between sites. The linear size is then equal 
to the total number of hollow and solid circles along L 
(see Fig. [T]) . At the i-th iteration, the linear size is given 
by Li = 2* + l for the arrowhead lattice and Li = 3 x 2 4_1 
for the crab lattice. Indeed, one observes that the di- 
mension of the lattice obeys the following relation with 
the total number of sites and linear size of the lattice: 
d g = lim^oo In N / In L t . 

We implement the Abelian Manna model on these two 
fractal lattices that have the same dimension but pos- 
sess different microscopic structures. Let us denote the 
variable Zj, which is a non-negative integer, to be the 
"height" or the number of particles at lattice site j. The 
lattice is first initialized with Zj = for all sites. The 
value of Zj is then evolved according to the following al- 
gorithm. When the system is in a stable configuration, 
i.e. Zj < 1 for all sites, the external drive is implemented 
by picking a site j at random and increasing its zj by 
one unit. At this juncture, an avalanche might occur in 
the following manner: as long as there exist any j with 
Zj exceeding the threshold z c = 1 (an "active" site), pick 
one of them at random, say k, and reduce Zk by two. At 
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FIG. 1: Sierpinski arrowhead (top) and crab (bottom) lattice 
at 4 th iteration with definition of linear size. In the figure, 
L = 17 for arrowhead and L — 24 for crab lattice. 

the same time, pick two of its nearest neighbours, say 
k' , randomly and independently, and increment zy by 
one unit. This procedure constitutes a toppling, which 
in the bulk is conservative, i.e. the total Zj remains 
unchanged by the bulk topplings. In one dimension, ev- 
ery bulk site has two neighbours. On the fractal lattices 
described above, a bulk site can have either two or three 
nearest neighbours. Note that particles can leave the 
system at the two end sites (labeled as 1 and 3* + 1 re- 
spectively). Owing to the Abelianess of the model, the 
statistics of the avalanche sizes is independent from the 
order of updates. This is different for time-dependent 
observables such as the avalanche duration, which is not 
studied in the following. 

An avalanche ceases as soon as Zj < 1 for all j (qui- 
escence). The size s of the avalanche is measured as 
the number of topplings performed between two quies- 
cent configurations. The probability density V (s) for an 
avalanche of size s to occur is expected to follow simple 
(finite size) scaling 

V{s) = as-^{^) (1) 

asymptotically in large s > so with lower cutoff so, lin- 
ear system size L, non-universal metric factors a and b 
and universal exponents r and D. The universal scaling 
function Q decays, for large arguments, faster than any 
power law, so that all moments 

(««) = JdsV{s)s q (2) 

exist for any finite system. Provided that q + 1 — t > 
one can easily show that 

(««) oc L D( - q+1 - T ~) . (3) 



The results presented in the following are based on 
Monte Carlo simulations for four different system sizes, 
corresponding to four different levels of iteration i — 
3,4,5,6, containing N = 28,82,244,730 sites and with 
linear size L = 9, 17, 33, 65 for the Sierpinski arrowhead 
lattice and L = 12,24,48,96 for the crab lattice. In all 
four cases, 10 s avalanches were triggered and the data 
recorded in regular intervals of 10 6 avalanches. Station- 
arity was verified by inspecting avalanche size moments 
and the transient was determined to be shorter than 5T0 4 
avalanches. Errors for the moments are derived from a 
Jackknife estimator [T^J [T3] of the variance based on the 
moments taken in each set of 10 6 measurements. The 
scaling exponents were determined by a nonlinear least 
square fitting |14j of the avalanche size moments against 
the linear size of the lattice: 

(s*) = (a, + \ + ||) L D ^^ . (4) 

The exponents derived in this procedure were r = 
1.170(5) and D = 2.792(2) for the arrowhead lattice; and 
t = 1.153(4) and D = 3.026(2) for the crab lattice. 

In addition, the scaling behaviour can be probed in a 
data collapse, as illustrated in Fig. [2] 
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FIG. 2: Data collapse of binned data for Sierpinski arrowhead 
lattice. Preliminary statistics of large system size (7 th itera- 
tion, L = 129) are also included to test the consistency of the 
data. The plots confirm the estimated values of the critical 
exponents. 

With the critical exponents determined, one can im- 
mediately verify the scaling law of avalanche size distri- 
bution. From hypercubic lattices it is well known [15] 
that the first moment of the avalanche size is given by 
the expected number of moves that a random walker per- 
forms on the given lattice before it reaches the boundary 
and leaves, i.e. by its residence time. This is essentially 
because of bulk conservation: In the stationary state 
one particle leaves the system for every particle added 
(avalanche attempt) and the average number of moves it 
performs during its residency is exactly twice the average 
number of topplings occuring in the system per particle 
added, which is the avalanche size. Regardless of the 
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specifics of the boundary, i. e. regardless of whether only 
two sites are dissipative or all sites along the perimeter, 
the first moment normally scales with the linear size of 
the lattices squared, D(2 — r) = 2, independent from the 
dimension of the hypercubic lattice. This is easily under- 
stood, as the time and thus the total number of moves 
performed by a random walker scales quadratically in the 
linear distance traversed. 

For the fractal lattices, it is obvious that D(2— r) is not 
equal to 2 since it is 2.317(8) for the arrowhead lattice 
and 2.564(6) for the crab lattice. We will now show that 
the scaling law has in fact changed to D(2 — r) = d w , 
where d w is the fractal dimension of random walk on the 
lattice. The scale law D(2 — t) = d w remains true for any 
lattice regardless of dimension and microscopic details. 

We will now calculate d w for the Sierpinski arrowhead 
lattice and the crab lattice by using the first-passage time 
method |16j . Due to the nearest neighbor structure of 
the fractal lattices, the calculation is not coarse-grained 
renormalization-like, but rather, it is carried out by con- 
sidering every single edge that connects between the two 
end sites on the lattice. 

We label the sites of the lattice sequentially with j = 

(i) 

1, 2, • • • , 3* + 1 and denote by Tj the average time for 
the random walker to exit through the end sites from site 
j (i is the number of iterations of the lattice) . Also, we 
denote by t the traverse time from one site to its nearest 
neighbor. 

Let us define E^ 1 ' to be the set of extended sites at the 
i-th iteration. On the lattice, each normal site j has two 
nearest neighbors j ' — 1 and j + 1. If j is an extended 
site, it has in addition a third nearest neighbor j*. If the 
walker is at site j, we have 2 situations: 

If j ^ the walker has only 2 options to choose 

from: go to j — 1 or go to j + 1. The probability for each 
choice is 1/2. We have 



-(t + T W 
2 V + 



t ± -T 



-T 



(5) 



If j G E^ , the walker has up to 3 options to choose 
from: go to j — 1, go to j + 1 or go to j* . The probability 
for each choice is 1/3. Thus we have 
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Since this gives a system of linear equations, we can write 
them in the form of a matrix equation: A^'T^' = B^\ 
where is the column vector of , is the degree- 
normalised adjacency matrix, while B^ is a column vec- 
tor contains t in every entry. As A^ is invertible, we can 
solve for T« = B®. By defining the new vari- 

able: — T^ %+r> /T^ , we obtain the results as shown 
in Table U for the arrowhead lattice and the crab lattice. 



TABLE I: Convergence of S[ as the number of iterations 



i increases. After i — 7, we have ~ 5 for the Sierpinski 
arrowhead lattice and ~ 6 for the crab lattice. Other 
Sj 8 are found to converge to the same limits. 
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Arrowhead 11 5.5076 4.7553 4.8806 4.9610 4.9893 4.9974 
Crab 10 7.1417 6.9528 6.3729 6.1014 6.0303 5.9727 
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FIG. 3: Convergence of as number of iterations i in- 
creases. After i = 7 we have S{" = 5 for arrowhead lattice 
and ~ a r>+v,«^ cf 1 ) 

same limits 



6 for crab lattice. Other S, • s converge to the 



We observe that for both lattices Sj converges to 
some value S* as i increases. For the arrowhead lat- 
tice, S* — 5, while for the crab lattice, S* w 6. On going 
from the i-th to (i + l)-th iteration, the lattice size in- 
creases by a factor Li + \/Li = 2, while the escape time 
increases by a factor T^ +1 ) /T« = S* for large i, so that 
d w = In Sj/ In 2. By comparing the calculated d w to the 
estimated value of D(2 — r) from simulation, we can see 
that they are in good agreement (see Table [IT| . 

One surprising conclusion from the above results is 
that for two lattices with the same dimension but dif- 
ferent microscopic structure, the critical exponents D 
and r can both be different, which suggests that on 
fractal lattices, the critical exponents depend not only 
on the dimension but also on the microscopic details 
of the lattice. In addition, our results have validated 
the scaling relation D(2 — r) = d w for two fractal lat- 
tices with different d w . Since on lattices with inte- 
ger dimensions d w = 2 regardless of the lattice's struc- 
ture, this scaling relation generalizes the standard ver- 
sion D(2 — t) = 2 known for the hypercubic lattices. 



TABLE II: The calculated d w and the estimated D{2 - 
for the Sierpinski arrowhead lattice and the crab lattice. 



Lattice Arrowhead 



Crab 



d w 2.322 2.578 

D(2-t) 2.317 ±0.008 2.564 ±0.006 
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TABLE III: Data for (d g ,d w , D) from 6 different lattices and 
the calculated Dtriai from Eq. Q. 

Lattices ID 2D 3D 4D Arrowhead Crab 

d g 1 2 3 4 1.58 1.58 

d w 2 2 2 2 2.32 2.58 

£> 2.2(1) 2.73(2) 3.36(1) 4 2.792(2) 3.026(2) 
Dtriai 2.10(7) 2.73(7) 3.36(7) 4.00(7) 2.78(7) 3.03(7) 



the calculated d w and the estimated D(2 — t) in Table [III 
They seem to be caused by finite size corrections, which 
are further suppressed at 7 iterations and above, as ob- 
served in preliminary data not included in the main anal- 
ysis above. In the presence of strong finite size correc- 
tions and high accuracy measurements of the moments, 
the estimates for D and r are sensitive to the choice of 
the fitting function, which is constrained by the number 
of data points (i.e. system sizes) available, but needs to 
contain as many correction terms as possible to account 
for the accurate data. Our choice Eq. Q reflects the de- 
sire to reduce the sensitivity of the estimate on the initial 
values. 

It is important to note that there is a level of ambiguity 
in the finite size scaling in fractal lattices, because due to 
its highly irregular nature, there is, a priori, no unique 
way of increasing the lattice size of a fractal [21] . At a 
given level of iterations in order to increase the lattice size 
further, one might either proceed by iterating the fractal, 
or use the given fractal to tessellate the hypercubic lattice 
of appropriate (embedding) dimension. One might argue 
that finite size scaling is of course sensitive to that choice 
and as a result, generates asymptotically either the ex- 
ponents of the fractal lattice or of the embedding space. 
However, in ordinary critical phenomena, there are cases 
[21] where the (effective) critical point and even the scal- 
ing functions change with the level of iteration i. In the 
current context that translates to, for example, the am- 
plitude A q = a q + b q /L + c q /L 2 in Eq. Q to acquire a 
dependence on i, which might distort the resulting esti- 
mates. The exponents derived above can thus be seen 
only as effective exponents of a fractal lattice. 



The last unexpected conclusion pertains to the relation 
between the dimensions d g and d w , and the critical ex- 
ponent D. We found that they obey the general linear 
relationship: D — ad g + /3d w + 7. This was uncovered by 
plotting the following six points in 3D parameter space: 
(d g ,d w ,D) = (1,2, 2.1), (2, 2, 2.73), (3, 2, 3.36), (4, 2, 4), 
(1.58, 2.32, 2.79), (1.58, 2.58, 3.03), which correspond to 
the linear chain [THj , the square pj |T7] , the cube [T51 [H] , 
the hypercube [50], the arrowhead and the crab lattices 
respectively. It is interesting that while the first four 
points due to the hypercubic lattices are found to lie on 
a straight line, all six points together make up a plane 
instead of a tetrahedron. We have determined the co- 
efficients a = 0.632 ± 0.003, (3 = 0.980 ± 0.014 and 
7 = —0.492 ± 0.034 from the six points, which leads to 
the following relationship: 

D ~ D trial (d g , d w ) = 0.632d 9 + 0.980d w - 0.492 . (7) 

A comparison of the value of D determined by Eq. Q 
and that from the simulation is shown in Table IIII1 
Finally, we comment on the slight mismatches between 
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